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Abstract 

Uncertainty in the determination of emission angles of signal and idler photons is investigated. The 
role of divergence of the pump beam, width of the pump spectrum, and size of the crystal is elucidated. 
Experimental data obtained by an iCCD camera are in a good agreement with a numerical model that 
provides angular and spectral characteristics of the signal and idler fields. 
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1 Introduction 

Light emitted from parametric fluorescence in a nonlinear crystal is composed of photon pairs. During 
the emission frequencies and emission directions of two photons comprising a pair are determined by 
the laws of energy and momentum conservation, as has been recognized already in year 1968 [1]. For 
this reason, there occurs a strong correlation (entanglement) between properties of the signal and idler 
photons. In an ideal case of infinitely long and wide nonlinear crystal and monochromatic plane-wave 
pumping, a plane-wave signal photon at frequency lo s has just one plane-wave idler photon at frequency 
U{. Emission angles of two photons are given by momentum conservation that forms phase-matching 
conditions. Possible signal (and similarly idler) emission directions lie on a cone which axis coincides 
with the pump-beam direction of propagation. In real experimental conditions, blurring of emission 
directions occurs because of crystals of finite dimensions [2,3], pump-beam divergence [4,5] as well as 
pulsed pumping [6,7]. Accepting the approximation based on a multidimensional gaussian two-photon 
wave- function an analytically-tractable theoretical model has been developed [8,9]. Recently, theoretical 
models have been generalized to photonic [11] and wave-guiding [12] structures. 

Here we extend the previous investigations of photon-pair properties by studying correlation areas of 
the signal and idler photons. Experimental results are compared with a theoretical model that considers 
Gaussian temporal spectrum and elliptical pump-beam profile. 

Model of parametric fluorescence is briefly described in Sec. 2. Details of the experimental setup are 
revealed in Sec. 3. The obtained experimental results are discussed in Sec. 4. Sec. 5 provides conclusions. 
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2 Model of parametric fluorescence 



The process of parametric fluorescence is described by the following interaction Hamiltonian H m t 
[2,10]: 

H int (t) = e o ydrx (2) :Ej +) (r,t)EH(r,t)Ei" ) (r,t)+H.c., (1) 
v 

where Ep + ^ is the positive-frequency part of the pump-field electric-field amplitude, whereas e! ^ (e| ) 
stands for the negative-frequency part of the signal- (idler-) field electric-field amplitude operator. Symbol 
X^ 2 ' denotes the second-order susceptibility tensor and : is shorthand for tensor reduction with respect 
to its 3 indices. Susceptibility of vacuum is denoted as so, interaction volume as V and H.c. substitutes 
Hermitian-conjugated terms. 

Parametric fluorescence of type-I in an LHO3 crystal oriented such that its optical axis was perpen- 
dicular to the z axis of fields' propagation was investigated. The pump field was polarized vertically 
(it propagated as an extraordinary wave) whereas the signal and idler fields propagated horizontally 
polarized (as ordinary waves). Under these conditions, a scalar theory is sufficient for the description. 
Because of low signal- and idler-field intensities, the Schrodinger equation can be solved only to the first 

order [7] and this solution then provides the fourth-order correlation function of the signal- and idler-field 

(2) 

electric-field amplitudes G s ■ in the form: 

G^(k s0 ,lqo)= f d 3 k s [ d 3 k i |d(k 5 )d(k i )| 2 |5(k a ,k i )| 2 J (2) 

where d(k s ) [d(ki)] is amplitude transmissivity of a detector for the signal [idler] field plane wave with 

(2) 

wave vector k s [kj]. The correlation function G s i gives the number of photon pairs generated with a 
signal-photon wave vector k s in the area defined by Ak s around k s o and the idler-twin wave vector kj in 
the area around kjo described by Akj. Correlation function S introduced in Eq. ([2]) can be expressed in 
the form 

S(k s ,ki) = C n j d 3 k p Ej+\k p )5(LU s + uj> - uj p ) f d 3 rexp[-i(k p -k s -ki)r] (3) 

v 

that reflects the role of the pump beam in the determination of the signal- and idler-field properties. 
Physically, it gives the probability amplitude of simultaneous emission of a signal photon with wave 
vector k s together with an idler photon having wave vector kj. Symbol C n in Eq. ([3]) stands for a 

normalization constant. Integrals occurring in the expressions in Eqs. ([2]) and (J3j) are complex in general 

(2) 

and require numerical approach. Program that determines the correlation function G s [ for any type-I 
parametric fluorescence has been developed to obtain results needed for comparison with the obtained 
experimental data. 

3 Experimental setup 

In our experiment, a negative uniaxial crystal made of LHO3 (manufactured by EKSMA Optics) 
was cut for non-critical phase matching (optical axis perpendicular to the pumping beam). Two crystal 
lengths were used L cr =2 mm or 5 mm both for cw and pulsed pumping. In cw case, a semiconductor 
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Figure 1: Experimental setup used for the measurement of angular uncertainties: (a) Whole setup using 
both pumping options and diagnostics of the pump beam (see the text for detailed description); (b) Part of 
the setup mapping propagation angles behind the nonlinear crystal to the positions on the photocathode. 

laser Cube 405 (Coherent) delivered 31.6 mW at 405 nm with spectral bandwidth AA =0.59 nm (all 
widths are given as half- widths at 1/e of the maximum) . For pulsed pumping the second-harmonic field 
of an amplified femtosecond Ti:sapphire system (Mira+RegA, Coherent) lasing at 800 nm was used. The 
pulses were ~250 fs long at the fundamental wavelength. At a repetition rate of 11 kHz the mean SHG 
power was 2.5 mW at the crystal input. Spectral bandwidth could be adjusted between 1.7 and 2.6 nm 
using fine tuning of the SHG process. The SHG beam was separated from the fundamental beam using 
a dispersion prism (see Fig. [[]) . 

Divergence of both pumping options was controlled using interchangeable converging lens LI or beam 
expander (BE2X, Thorlabs). Focal lengths fn in the range from 30 to 75 cm were used. Distance 
between the lens LI and the nonlinear crystal was set such that the beam waist was placed behind the 
crystal, Zli < fn. Spatial spectrum of the pump beam was measured by a CCD camera (Li085M, 
Lumenera) placed at a focal plane of a converging lens L3. From the spatial spectra, vertical and 
horizontal projections of the beam waist Wqi, I = x,y, were determined. Temporal spectrum of the 
pump beam was measured by a fiber spectrometer (HR4000CG-UV-NIR, Ocean Optics) placed behind 
the crystal. 

In our experiment we focused on degenerate photon pairs A s o = = 800 nm that were emitted on 
a cone layer with apex angle 33.4 deg behind the crystal. One section of the cone layer was captured 
directly by the detector, the opposite section was directed to the detector using a high-reflectivity mirror 
(see Fig. [2h,). The detector was composed of an iCCD camera with image intensifier (PI-MAX:512-HQ, 
Princeton Instruments) preceded by a converging lens L2, one narrow-bandwidth and two high-pass edge 
filters. The purpose of lens L2 is to map the photon propagation direction angles to positions at the 
photocathode (see Fig. EJa). Several lenses L2 with focal lengths Jl2= 12.5, 15, or 25 cm were used in 
different variants of the experiment. Lens L2 was placed such that the photocathode lied in its focal plane. 
Filter 11 nm wide was centered at 800 nm. Edge filters (Andover, ANDV7862) blocked wavelengths below 
666 nm while their transmittance at 800 nm reached 98%. 

Photosensitive area of the detector in the form of a rectangular shape 12.36 mm wide was divided 
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Figure 2: (a) Cumulative registration of photons produced by 20,000 consecutive pump pulses in one 
camera frame. Signal and idler strips image small sections of the cone layer (note slight curvature of the 
strips, the curvature is oriented in the same sense for both strips because the idler strip is reflected from 
a mirror), (b) A typical frame containing photon detection events caused by one pump pulse. 

into 512 x 512 pixels. Resolution of the camera was limited to 38 ^m (FWHM) mainly due to imperfect 
contrast transfer in the image intensifier. To speed-up the data collection we decreased the resolution 
even further by grouping 4 x 4 or 8 x 8 pixels into one super-pixel in the hardware of the camera. Typically 
we captured several tens of camera frames per second. Quantum efficiency of the detector was estimated 
to 7% including all the components between the nonlinear crystal and photocathode. 

To measure the signal-idler correlation function given in Eq. ([2]) we defined three regions-of-interest at 
the iCCD photocathode (see Fig. [2|). Two of them capture signal and idler photons and their widths are 
given by bandwidth filter and lens L2 focal length. The third strip serves for monitoring the noise level. 
In a pulsed regime a 10 ns long gate of the camera was used synchronously with laser pulses. Considering 
cw pumping the camera was triggered internally with a gate lasting 2 /j,s. Pump-field intensity was set 
so that the average number of photon detection events per signal/idler strip was much lower than the 
number of super-pixels in each strip (~ 5000). This made the probability of detecting two photons in a 
single super-pixel negligible. 

Not all detection events were due to signal or idler photons. We made a detailed analysis of the noise 
and found out that 1.82% of detections came from noise the majority of which (9/10) were red photons 
coming from fluorescence in the crystal. The rest were scattered pumping photons and only 1.6% of the 
noise was caused by dark counts. 

Evaluation of the detection frames allowed us to compute the experimental correlation function in 
the horizontal (X) and vertical (Y) coordinates at the camera (see Fig. [2b) as follows: 

N M p L p 

g (X S , = EWT 6 W - X s) S (Xf ~ Xj) , (4) 

p=l m =l 1=1 p p 

where p is an index of the frame (N gives the number of frames) and m (I) counts signal (idler) detection 
events [up to M p (L p ) in the p-th. frame]. Symbol Xg l [Xf] denotes horizontal position of the l-th. 
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Figure 3: Results of a typical measurement sequence with pulsed pumping containing 327,600 frames, 
L cr = 5 mm. (a) Spatial spectrum of the pump beam as measured in lens L3 focal plane, (b) Intensity 
spectrum of the pump beam; diamonds are experimental values from spectrometer, solid line is a multi- 
peak Gaussian fit. (c,e) Experimental correlation functions g(Xs,Xi) and g(Y~s,Yj) evaluated according 
to Eq. (j4|). (d,f) Cross-sections of functions g(X$,Xj) and g(Yg,Yi) along the signal axis averaged over 
all registered idler values (open rectangles) and the corresponding theoretical curves G given by Eq. (J2]) 
obtained from the numerical model (solid lines). These plots are cast in angular units j s ~ X s //l2, 
(is « V, ./,.,. • • - 

detection in the signal [idler] strip in the p-ih frame. Correlations in the vertical direction are determined 
similarly. In this way, all possible combinations of pairwise detection events are taken into account. 

4 Results 

A typical result for pulsed pumping based on a measurement sequence lasting three hours is given in 
Fig. 02 Spatial spectrum of the pump beam as measured in the focal plane of lens L3 is shown in Fig. [3^ 
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Figure 4: Dependence of vertical angular uncertainty (A/3 S ) on the bandwidth of vertical pump-beam 
waist W y o for pulsed (crosses and triangles) and cw (diamonds) pumping. Two different crystal lengths 
were used: 5 mm (crosses) or 2 mm (triangles and diamonds). Solid line is provided by the numerical 
model. Dash-dot horizontal line indicates the resolution limit given by camera super-pixel size. 



and the corresponding intensity spectral profile in Fig. [3b. Plots of the correlation functions g(Xs,Xj) 
and g(Ys,Yj) evaluated from 327,000 registered frames according to Eq. (jj]) are given in Figs. [3h and 
[3^. Pairwise correlated character of the fields leads to the diagonal patterns going from upper-left to 
lower-right corners of the plots. The diagonals have finite widths originating in spreading of signal 
propagation directions related to one fixed idler direction. We denote half-widths of the diagonals as 
AX s /2 and AY s /2 or, in angular quantities, as Aj s /2 ps AX s /(2f L2 ) and A(3 s /2 w AY s /(2f L2 ). Since 
these widths apparently do not change significantly in the X and Y ranges used for plotting Figs. Eh and 
[3k, we can increase the precision in determining A7 S and A/3 S by averaging over the range of measured 
idler coordinates. The averaged cross-sections along signal coordinates 7 S and j3 s are indicated by angle 
brackets () and are plotted in Figs. [3U andEF (open rectangles). Here, solid lines show the corresponding 
curves giving G(k s ,kjo) that are obtained from the numerical model using parameters of the pump 
beam as were derived from the measurements reported in Figs. [3ti and [3b- Excellent agreement of the 
model with experimental data is evident even though the pump-field spatial shape and its spectrum are 
nontrivial. 

We have performed systematic investigation of the dependencies of angular uncertainty (A/3 S ) on 
various parameters of the pump beam and crystal length. The role of vertical dimension of the pump- 
beam waist Wyo in the behavior of the uncertainty (Aj3 s ) is revealed in Fig. U] where values for both 
pumping options and two different crystal lengths are plotted on the top of the solid line given by the 
numerical model. We note that in this case and according to the model values of the uncertainty (A/3 S ) 
do not change with pump spectral width and horizontal beam- waist size W x o- We can see in Fig. H] 
that experimental points follow the theoretical curve except for the area of narrow diagonals where the 
experimental resolution was limited by the size of super-pixels. 
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5 Conclusions 



Thanks to the ability of the intensified CCD camera to register single-photon detection events and 
positions of their occurrence, we were able to measure angular uncertainties of the far-field signal-idler 
correlations in the process of parametric fluorescence in an LIIO3 nonlinear crystal. In particular, we 
have measured the dependence of vertical angular uncertainty on the vertical pump-beam waist size using 
both cw and pulsed pumping and two different crystal lengths. We have reached a very good agreement 
with the numerical model of parametric fluorescence that simulates the process for any polychromatic 
elliptical pump beam. 
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